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Kontsevich integral for knots and Vassiliev invariants 

P.Dunin-Barkowski* A.Sleptsov| A.Smirnov^ 



ABSTRACT 

We review quantum field theory approach to the knot theory. Using holomorphic gauge we obtain the Kontsevich integral. It is explained 
how to calculate Vassiliev invariants and coefliicients in Kontsevich integral in a combinatorial way which can be programmed on a computer. We 
discuss experimental results and temporal gauge considerations which lead to representation of Vassiliev invariants in terms of arrow diagrams. 
Explicit examples and computational results are presented. 
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1 Introduction 



It is believed that the path-integral representation for knot invariants arising from topological quantum field theory (QFT) 
gives the most profound and general description of knot invariants. Ideologically, it means that all possible descriptions of 
knot invariants can be derived from this representation by utilizing different methods of path-integral calculus. For example, 
the usage of certain non-perturbative methods leads to well known description of polynomial knot invariants through the 
"skein relations" jT]. The perturbative computations naturally lead to the numerical Vassiliev Invariants j2j. In the last case 
we obtain the formulae for Vassiliev invariants in the form of "Fcynman integrals" . For a recent comprehensive treatise on 
Vassiliev invariants see [3]. 

Despite beauty and simplicity of this picture many problems in the theory of knot invariants remain unsolved. Currently 
more mathematical descriptions of knot invariants are developed than can be derived from path-integral. One such problem is 
the derivation of quantum group invariants from the path-integral representation. The main ingredient in the theory of these 
invariants is the universal quantum i?-matrix defined for integrable quantum deformation of a Lie group. The appearance 
of the quantum groups in the path-integral representation looks mysterious and we lack the derivation of the corresponding 
i?-matrix (object with noncommutative matrix elements should appear from classical integral) from path-integral. These 
problems were discussed in details in [4, 5 . 

The second interesting problem (which as we believe is closely connected to the first one) concerns the combinatorial 
description of the numerical Vassiliev invariants of knots. At the moment, there exist three different description of Vassiliev 
invariants: through the generalized Gauss integrals [6j, through Kontsevich integral [7], and finally there are combinatorial 
formulae for Vassiliev invariants of orders 2,3 and 4, [8]- [11]. The first two descriptions can be easily derived from the 
path-integral but (surprisingly) we lack such a derivation for combinatorial formulae. The Gauss integral representation for 
Vassiliev invariants comes from the perturbative computations of path-integral in the covariant Lorentz gauge. Similarly, 
usage of non-covariant holomorphic gauge (which is sometimes referred to as light-cone gauge) computation of path-integral 
leads to Kontsevich integral [12\ . The aim of this paper is to show that the universal combinatorial formulae for Vassiliev 
invariants and quantum group invariants can be derived from path-integral in the temporal gauge. 

The main difference between the non-covariant gauges (temporal or holomorphic gauge) and covariant Lorentz gauge is 
that the Feynman integrals in there gauges can be naturally "localized" , in the sense that only a finite number of special points 
on the knots contribute to the integrals. In the Lorentz gauge the Feynman integrals have a form of multiple 3-d integrals (see 
([6|,([7]) as an example) and all points of the knot enter the integral equally. On the contrary, in the holomorphic gauge only 
special points contribute to the Feynman integrals. Summed in all orders of perturbation theory these contributions lead to 
definition of polynomial knot invariants through simple crossing operator and Drinfeld associator with rational zcta-function 
coefficients [13] . We discuss this localization process in detail in Section [s] 

Before we pass to temporal gauge consideration let us briefly describe Section |4j where the computation of Vassiliev 
invariants is discussed. These results are used in Section [5] for numerical studies. To compute Vassiliev invariants we use 
technique developed in previous Section [Sj With the help of it we can calculate coefficients in Kontsevich integral expansion 
in chord diagrams. However, these coefficients are not invariants due to chord diagrams being not independent. We choose a 
basis among the independent ones with the help of so-called trivalent diagrams. Thus we can find Vassiliev invariants. In the 
rest of the Section we describe two polynomials in Vassiliev invariants up to level 6 vanishing on a large class of considerable 
knots. 

The temporal gauge is distinguished among all gauges. The Feynman integrals here have ultra-local form - only crossing 
points of two-dimensional projection of the knot contribute to the answer. This is exactly what happens in the quantum group 
description of knot invariants where the crossing points contribute as a universal quantum i?-matrix. On the other hand, 
the combinatorial formulae for Vassiliev invariants also are based on the information from these crossing points only. All 
these observations make it natural to argue that the Feynman integrals arising from the general path-integral representation 
in the temporal gauge should give the universal combinatorial formulae for Vassiliev invariants, and summed in all orders of 
perturbation theory they should give the perturbative expansion of the universal quantum i?-matrix for /i-deformed gauge 
group of the path integral. These questions and structures arising from path-integral representation of knot invariants in the 
temporal gauge are discussed in Section [5] 

As a side remark, recently new polynomial invariants of knots appeared, the so-called "superpolynomials" [131 [H], which 
are generalizations of HOMFLY polynomials. However, there is still no related QFT-like formulation as the one discussed in 
the present paper. Probably this interpritation should involve something like "Chern-Simons theory with a quantum group" 
which is still not constructed (for an attempt see [16]). 
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2 Knot invariants from C hern- Simons theory 

The incorporation of Vassiliev invariants in the path-integral representation is clear from the following picture. Let A be a 
connection on taking values in some representation i? of a Lie algebra g, i.e., in components: 

A = A^{x)T'' dx\ 

where T'^ are the generators of g. Let curve C in give a particular realization of knot K . Consider the holonomy of A 
along C, it is given by the ordered exponent: 

r(C,A) = Pexp j)A^l + ^A°(x)T" + j A'llixi) j ^^^^(xz) r°^T°^ + ... 

C C CO 

The Wilson loop along C is a function depending on C and A defined as a trace of holonomy: 

WR{C,A)=tvRV{C,A) 

According to [T] there exists a functional Scs{A) (we write it down explicitly later) such that the integral averaging of the 
Wilson loop with the weight exp ^ — '^S{A)^ has the following remarkable property: 



'Wr{K)) = - / i^Aexp 



{A))Wr{C,A) 



(1) 



where 



J DA exp(-^5cs(v4) 



i.e. the averaging of Wr{C,A) with the weight exp ^ — ^S{A)j does not depend on the realization C of the knot in M'^ 

but only on the topological class of equivalence of knot K (in what follows we will denote the averaging of quantity Q with 
this weight by (Q)) and therefore, {W{K) ) defines a knot invariant. 

The distinguished Chern-Simons action giving the invariant average ([T]) has the following form: 

(2) 



ScsiA)= I tiiAAdA+ -AaAaA) 



If we normalize the algebra generators as tr(r"r'') = (5°^ and define the structure constants / of algebra g as [T'^jT^] = 
fabc T'^ then the action takes the form: 



Formula ([I]) is precisely the path integral representation of knot invariants. It is believed that all invariants of knots can 
be derived from this expression. Let us outline the appearance of Vassiliev invariants in this scheme. Obviously the mean 
value ( W{K) ) has the following structure: 

{W{C,A))^{Y, (fdxi fdX2... f dxnA'''{xi)A'''{x2)...A''Hxn) tr{T'''T''\ . .T'^-)) = 

n=0 J J J 

oo „ „ „ oo Nn 

= Y,fdxi / dx2... / dxn{A-Hxi)A-Hx2)...A-'ixn)) tr{T-'T-\..T''-) = Y.J2 ^n,™ G„,„ (3) 



ri=0 



n— m— 1 



From this expansion we see that the information about the knot and the gauge group enter < W{K) > separately. The 
information about the embedding of knot into M.^ is encoded in the integrals of the form: 



Vn 



dxi / dx2... / dXn{A'''{xi)A''^{x2)...A''HXn)) 
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and the information about the gauge group and representation enter the answer as the " group factors" : 

Gn,m ~ tr(r''iT°^..r''") 

Gk,m are the group factors called chord diagrams with n chords. Chord diagrams with n chords form a vector space of 
dimension A^„. Despite {W{K)) being a knot invariant, the numbers Vn^m are not invariants. This is because the group 
elements Gn.m are not independent, and the coefficients Vn^m are invariants only up to relations among Gn,m- In order to 
pass to the the knot invariants we need to choose in the space of group elements some basis of independent group elements 
and expand the sum ^ in this basis. Finally as it was proven in [IJI we arrive to the following infinite product formula: 



: W{K) )^dR U n exp(fi"V„.^g„,™) 

n— 1 m— 1 



(4) 



where Af^ is the number of independent group elements of degree n and V„^„i are the coefhcients of independent group factors 
Qn.m which are nothing but the Vassiliev invariants of degree n; dn stands for the dimension of representation R. Note that 
the number of independent Vassiliev invariants of degree n is given by the number of independent group elements Afn- We 
list the first several values of Afn ■ 

n 1 9 S 4 .f^ fi 

(5) 



n 


1 


2 


3 


4 


5 


6 




1 


1 


1 


2 


3 
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Formula Q along with table ^ means that the expansion of {W{K) ) up to order 3 is the following one: 



:w{K): 



1 + Wi, 1^1.1 + h' 



3! 



"^"1 1^1,1 + Vl,lV2, 1^1, 1^/2,1 + V3,l^?3,l 



Note, that here the relations between group factors are taken into account. 

QFT provides several techniques for computing coefficients of /i-expansion for < W{K) >, each technique leads to some 
formulae for Vassiliev invariants. The straightforward way is to use the perturbation theory for covariant Lorentz gauge 
diAi = 0. Standard quantization technique for Lorentz gauge considered in 6 leads to the following Feynman integral 
formulae for the first two Vassiliev invariants (integrals are taken along the curve representing the knot): 



dx\ 



^'^^'^^■'^Ixi-x^P 



(6) 



V2,i ^^1 J dxi. j dxl J dx'^epjgEkis- 



xa - X2\^ \X3 - XiP 



c 

-^-Idx' 7dx' Idx'^ /dx^eP-e e- , (^4 - ^i)" (2:4 - 2:2)" (^4 - 0^3)* 
^fdx.J dx, J dx, J dx^e e^p^e.rnekst _ _ _ (0 

C R3 

The first integral here is the so called Gauss integral for self linking number. In general, perturbation theory in Lorentz gauge 
provides the theory of Vassiliev invariants in terms of rather sophisticated 3 — d " generalized Gauss integrals" . Therefore we 
will consider different choices of gauge in the following Sections. 



3 Holomorphic gauge 

The aim of this Section is to show that the universal Kontsevich integral for knots can be obtained from the general path 
integral representation. In fact, as it was explained in [T^ the Kontsevich integral appears as average value ( W{K) ) computed 
in the so called holomorphic gauge. To show it, let us start with decomposition of three-dimensional space = M x C, i.e. 
we pass from the coordinates {xo,xi,X2) to {t,z,z) defined as: 

t — Xq, z — Xi + iX2, z = Xi — iX2 
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then the for differentials and dual bases we have: 



dz = dxi + idx2, dz = '^i'^^i ~ ^'^^2) 



dz ~ dxi — idx2, dz 



The gauge field takes the form: 



where: 



A-^dx' = A^dx° + ^A^^dz + ^A^^dz 



Aa _ Aa _ ■ Aa Aa _ Aa , ■ Aa 
■^z — ^1 '^2 J ^2 — ^1 + '^2 

As usual, the path integral ([T]) is not well defined due to gauge symmetry, to make it convergent we need to fix an appropriate 
gauge condition. The holomorphic gauge is defined by the following non-covariant condition 



A-z^Q 



(8) 



The main feature of this gauge is that the initially complex, cubical Chern-Simons action becomes pure quadratic and the 
resulting path integral becomes Gaussian. This allows one to compute it perturbatively utilizing the Wick theorem. Indeed, 
in the gauge ([s]) the cubical part of the action ([2| vanishes: 



A^A^A\ 



and we end up with the following quadratic actior[^ 




(9) 



where e™" is antisymmetric, m,n € {t,z} and e*^ = 1. The main ingredient of perturbation theory with action ([9| is the 
gauge propagator defined as the inverse of quadratic operator of the action: 



[A'^^ih, Zi, Zi) Ai{t2, Z2, Z2)) 



gab 



h &{ti-t2] 

27ri Zi — Z2 



(10) 



To find it we need the following simple fact about the operator dz'. 



27ri z 

To prove it, we note that the inverse of Laplace operator A = dzdg (its Green function) on the complex plane is given by 
the logarithm functiorr] 

1 111 

{dzdzy^ = —\og{zz), therefore d^^ = —dzlog{zz) ^ — - 
Zm Zm IiTi z 

Let us show that average value {W{K)) in the holomorphic gauge coincides with the Kontsevich integral for the knot. To 

proceed we introduce an orientation on a knot in M^" . Obviously, if we pick some point p on the knot then K\p is topologically 

a segment Ip = (0, 1). Orientation on Ip naturally defines orientation on K\p. We will denote this orientation using symbol 

o, for example we can compare two points xi and X2 as o{xi) < 0(0:2). The Wilson loop is given by ordered exponent which 

has the following form: 

/CO n 
Al{x)T''dx^'^ =^ / \{ dxi'' Al\{xi)All{x2)...A';;^S^^)iY (t''^T'\..T'' 

o(x\)<o\xi)<...o(Xn) 



^Of course, the careful gauge fixing involves the Faddeev-Popov procedure which leads to an additional ghost term in the action. Fortunately, 
in the holomorphic gauge the ghost fields are not coupled to the gauge ones and therefore can be simply integrated away from the path integral. 
In this way we again arrive to the quadratic action (|9|. 



■^Note that 9^ ^ is defined up to any holomorphic function, as they are in the kernel of dz 



2m z 



The trick with the Laplacian consists of the following: we restrict the operators to the space of functions with absolute values sufficiently fast 
decreasing at the infinity. The only holomorpfic function with this property is f{z) = 0. 
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For the average ( W{K) ) we get: 

oo 



/n 
I — 1 



" ^ o(xi)<o{x2)<...o(Xn) ^ ^ 



As the action of the theory is quadratic, the average of n fields is not equal to zero only for even n, moreover the Wick 
theorem gives: 



E 



(11) 



)fej2)...(vj,.) e^': 



fc=l 



Here the sum runs over the set of all pairings P2n of 2n numbers. An element of this set has the formp = ^(ii, ji)(«2, j2)---(*rt, Jn)^ 
where ik < jk and the numbers ik,jk are all different numbers from the set {1, 2, 2n}. If p e is a pairing then we 
define a function Up on a set {1, 2, 2n} as follows: 

(^piik) = k, (^pijk) = ik 



i.e. it returns a minimum number from the pair {ik,jk)- The Wick theorem (11) gives: 

oo „ n 

^ „ ^ D. I 1 



n—0 I \^ f \^ / \ p£P2n k—1 



From the propagator (10) we have: 



dx^:" dx'/^" {x.^^^ (x.J) = ^ <5''-''- idz^.dt,, dz.jUj ^^h^^-hA 



therefore, 



oo .„ . n ^ 



where the pairing p = {^{ii, ji){i2, 32) ■■■{in, jn)^ and the group element Gp = tr |^T''''p(i)T''''p(^' ...T°''p'^"'y To obtain the 

Kontsevich integral we need to make the last step: integrate out the delta functions 6{ti^, —tj^,) from the integral. To do it, it 
is convenient to parametrize the knot by the "height" t. Such parametrization of a knot is regular at non-critical points with 
respect to t direction. In the interval between the critical points we have a well defined function z{t) which parametrizes the 
knot in M X C. Using this parametrization we rewrite the integral in the form: 



(27ri)" 

0(Zi)<0(22)<---0(2n) ^ 



k=l 



\ dti 



dti 



Here we should be careful not to forget the sign factor (—1)^^ where is the number of down-oriented segments between 
critical points on the knot entering the integral. On these segments the "height" parameter and the orientation of the knot 
are opposite, and therefore we should change dti — > —dti which finally results in factor (— 1)^^. Integrating the tj^ variables 
we obtain: 



00 n 

^w) = E7^ / E (-1)^^11'^*- 

o{zi)<o{z2)<...o{z„) 



dz,^{t,^) dZjJt,^) 



dt, 



dti 



1 



and finally, taking into account that dzji_{ti^) / dtj^, — dzjf,{ti^)/dti^ we arrive to the following expression: 



:w{K))= E 



{^niy 



I 



J2 (-ir A ^^^^^^^G, 



0(Z1)<0(Z2)<---0(Z„) 



k=l 



(12) 



The last expression ( 12 ) is the celebrated Kontsevich integral for knot K. 
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3.1 Localization of Kontsevich Integral 



3.1.1 Multiplicativity and braid representation 



The coefficients of Kontsevich integral ( 12 ) are given in terms of rather sophisticated meromorphic integrals. In this Section 



we describe the localization technique for computing the Kontsevich integral (KI) which gives a simpler combinatorial 



description of ( 12 ) [13j . The idea behind the localization is in fact very simple and based on miltiplicativity and factorization 



of KI and uses braid representation of knots. Let us discuss these properties separately. 




A 







: 


I" 




t2 






F mi 


m-2 



The multiplicativity of KI means that we can cut the knot in the finite 
number of parts and compute KI for these parts separately, then KI for the 
whole knot can be computed as an appropriate product of these separate inte- 
grals. To illustrate the idea, let us consider knot 3i (see [TH] for knot naming 
conventions) cut in three parts as in Figure [l] The Kontsevich integral com- 
puted for separate parts is not a number anymore but an operator of the form 
yf^Nin _^ y®No^t -^^here and No^t are the numbers of incoming and outgo- 
ing lines correspondingly. For example let us cut the knot in three slices (f4, fa), 
(^3,^2) and {t2,ti). Then the corresponding Kontsevich integrals are given by 



B 



operators A^-^^ 



where we write their indices explicitly 




to emphasize that they are finite-dimensional tensors. Lower indices correspond 
to incoming lines and the upper to the outgoing ones. The value of KI for the 
whole knot is the product of these tensors: 



E 



j\ili2 r)jlj2kik2 ^mim2 
jlj2 iii2mim2 feifc2 



Figure 1: Slices 



■ For the entire knot KI is a number and not a tensor because the knot is closed 
z 

and does not have any incoming or outgoing lines. Integral naturally brings us 
to the main idea of localization: if we are able to represent the knots as a union 
of finite number of some special "fundamental" parts, then we have to compute 
KI for these parts only. In order to compute the KI, one needs to choose these special parts in a way such that the KI 
for them would have the most simple form. Let us consider the middle slice (^2,^3) in Figure [T] Kontsevich integral B 
for this part is obviously equal to identity, this is because this piece consists of four lines parallel to t axis and the form 
(dzi ~ dzj)/{zi — Zj) in (12) vanishes as dz = Q on each vertical line: 



B = 1 



(»4 



Tjn\n2minn2 
kik2PiP2 



eni rn2 rmi rm2 
°ki°k2°Pl "P2 





Therefore, it is quite natural to use representations of knots " maximally extended" along the t direction, such that the KI 
takes the simplest possible form. Such a representation is the well known braid representation of knots: 

For any knot K there is a number n (not unique) such that it can be 
represented as a closure of some element ( not unique ) from braid group Bn ■ 

The meaning of this theorem is clear from the examples in Figure [2] Here 
the knots 4i and 82 are represented as closures of braids §291^9291 G -83 
and 929192 ^ 9i € correspondingly. The closure is the operation that 
connects the top of the braid with its bottom stringwise. There exists a 
simple combinatorial algorithm to construct braid representation h for any 
knot [IB]. 

The last thing we need is the factorization of Kontsevich integral. 
Let us introduce the distance between the strings in the braid. Then the 
following holds: if one can arrange the strings in the braid in a way that the 

distance between two groups of them is of order of different powers of some ^. „ -r^ . , . r , " 

. .1 i • . i_ -1 X .1 T^T rigure 2: Braid representation tor knots 

small parameter e, then these two groups give separate contributions to the KI. 

More precisely: Consider a braiding 6„ = bk ® 

bn-k G Bn- Let us assume that the sizes of bk and 

bn-k are much less (0(e)) than the distance between 

them. Then 

Klibn) - KI{bk) KI{bn-k) + 0{e) 

If bn — bk ^ bn-k then the first k strings of the braid 
do not cross the rest n — k strings. For example in 




Figure 3: Factorization 



Figure [3] the first three strings and the last two form 
two separate braids. If widths of these braids are ( 
^ e) much less the distance between them 1), then 

the factorization theorem implies that the KI is the tensor product of two KFs for each separate braid. 
3.1.2 Choice of associators placement 

Now, with the help of these properties we are ready to describe the representation of a given knot for which KI takes the 
simplest form. Let knot K be represented by closure of some braid 6 G S„, then we arrange the strings of the braid such 
that the distance between fc-th and k + 1-th strings is given by e~'^ where e is some small formal parameter, Figure H 
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2 


k - 1 


k 


k + : 




n 






























n 































-k-\-l -h -fc-1 Z 



Figure 4: R-matrix with associators 



i.e. the distance between adjacent strings increases with the number of string. Suppose, that at 
some slice (^4,^1) the braid b has the crossing of fc-th and fc -I- 1-th strings (Figure |4]). This crossing can 
be represented as a result of taking three consecutive steps: 

1. In slice (^4,^3) the fc-th string goes closer to the string fc -f 1, at the distance er'^^^. 

2. In slice (t3,t2) two strings are crossed. 

3. In slice (^2,^1) the fc-th string goes to its initial position at the distance e~'^ from fc -f 1 string. 

If we arrange the strings in this way, then the distance between fc-th and k + 1-th strings at the interval Figure 5 
(^37^2) is much less then the distance to any other string. Therefore, the factorization for KI gives the matrix 
following contribution for slice (i3,t2): 




Hum) - ^k,k+i - 

where R^^ is the so-called R-Matrix, i.e. KI computed for the neighborhood of the crossing point which corresponds to the 
configuration of strings in Figurejsj the sign depends on the orientation of the crossing. Similarly, at slices (t^, t^) and (i2, ^i) 
the distance between first and fc-th string is much less the the distance to the other n — k strings. Therefore, the factorization 
for the slice (i4,t3) gives: 



±1 



.1 



K 



±1 , 



1 



k+2--- 



In 



KL, 



where is the so-called associator, given by KI for the configuration of strings represented in Figure [Hj In summary, we 

see that the crossing of fc-th and fc + 1-th strings in braid representation of the knot gives the contribution to the KI of the 
form of operator R^f._^_i conjugated by operator ^'i,2,...,fc-i-i: 



= $ 



-k-l 



1 



®n—k—l 




Denote this tensor by Xk'- 



Xk — 5'l,2,...,/c-|-l Rk,k+1 



,fc+i- 



(13) 
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3.1.3 Formulas for R-Matrices and associators 



Explicit expressions for operators R and 5* can be easily derived from com- 
putations of KI for the corresponding pieces of the braid representation of a 
knot. The KI for neighborhood of a crossing point gives: 



R = exp ( fl 



(14) 



The calculation of $3 is explained in details in [T^, to describe it let us 
define: 

Q12 = (Xi 1, 5123 = 1 <8) fl 

then we have: 



$3 



1 



(813 



E 



k=2 



(—Y 



E E (-l)l'ilr(pi,gi,...p„,,9„) X 



Figure 6: Associator 



1>0 P>Oq>0 
|p| + |q| = fe 

'(p)='(q)=™ 



E (-1)'^' n 



!(r) = I(s) = n. 
0<r<p, 0<s^q 



^23 ^12 ^23 ■•■^12 ^23 ^12 (1^) 



kl 



Where p = (pi,P2, ■••,Pm) is a vector with positive integer components. The length of the vector ^(p) ~ m and |p| = YlPi- 
When we write p > q, it is understood as Pi > Qi for all i, and p > means that Pi > for all i. The coefficients 
t(pi, gi, ...prm Qm) SLve expressed through multiple zeta functions as follows: 



T{pi,qi, ...p,n,qm) = C(1^3'^l ^ l,l_^-^,'72 + 1, ■■■,qn + 1) 

pi-l P2-1 

such that for example t(1, 2) = C(3) and t(2. 1) ({I, 2). The multiple zeta functions are defined as: 

C(TOi,m2,...,m„) = ^ fcr""^2 

0<fci<fe2<. ..</£„ 



12 ...t-n 



Note, that t(pi, qi, ...p^, g™) = T{qm,Pm, ■■■qi,Pi) so that, e.g. C(lj2) = (^(2) and C(l:3) = C(4)- To define the associators 
$„ for n > 3 we need the coproduct operator A : U{g) — > U®'^{g). Its action on the generators of the universal enveloping 
algebra U{g) is defined as follows: 



A(T") = l(g)T'^ + T°«)l 

A{T''T'') = 1 r"T^ + T*- + T'' (g) + T°T^ 1 



(16) 



With the help of this operator the higher associators are expressed through $3 with the help of the following recursive 
formula: 



(17) 



Formula (17) basically means that, e.g. to obtain $4 from $3 one needs to make the following substitutions in the formula 
for $3: 

l^T^f^T" — ^ 1 ® 1 (g) «) T", 

(g) 1 (g) T" — > 1 (g) T" (g) T" (g) + (g) 1 (g) T" (E) T'' , 
rpaj^b (gjT" (gT'' — >l(g T'T'' (g (g T'' + T" (g T'' (g T'^ (g T'' + T'' (g (g T'' (g T'' + T'^T^ g) 1 (g T° g) 



and so on. That is, one has to symmetrize the first tensor component of each term of $3 over the first and the second tensor 
components of $4. 



One can derive (17 1 directly from Kontsevich integral consequentially for n = 3,4,..., but it is simpler to note that in 



Figure |4] the first k — 1 strings should give equivalent contribution to associator ^k+i because the distance between first and 
k— 1-th strings is much less then the "width of associator" (the distance between the first and the fc + l-th string). Therefore, 



operator A and the recursive procedure (17) have a clear physical meaning of "symmetrization" of contribution of the first 
A; — 1 strings to associator 
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3.1.4 Caps 

To complete our consideration we should also find contributions to KI coming 
from the bottom and the top of the braid's closure. Fortunately, we do 
not need new operators as these contributions can be represented through 
already introduced operators VPfc. 

Indeed, to make a closure of braid b E Bn we have to add to the braid 
n straight strings, such that the total number of strings is 2n and then 
connect {n — fc)-th string with the n + A:-th one for all k at the top and the 
bottom, for example as in Figure [7] In what follows we again imply that the 
distance between our 2n strings increases with the number of string, such 
that the distance between fc-th and k + 1-th ones is of order . Consider 
the following procedure, first the n-th string goes closer to n + 1 up to the 
distance e^" (slice (ii,i2) in Figure[7|) then we connect them by a "hat" of 



width 



The factorization theorem gives the contribution ^! ^ 



-1 

2. 



,n+l 



for 



this slice (due to the hat contribution being trivial). In the slice (^2,^3) we 
have 2n — 2 strings and we can iterate the procedure one more time which 
will give the contribution of the form ^+2- Finally, the whole top 

closure gives: 



-1 



-1 



-1 

l,2,2n-l * l,2,...n-l,n+2-- - ■^l,2,...ri+l 



(18) 




1 n - 2 n- — 1 n. fi + 1 n. -|- 2 n. -|- 3 2n. 



Figure 7: Cap 



It is assumed here that the indices corresponding to strings which terminate in a cap are contracted. 
Analogously for the bottom of the closure: 



l,2,2n-l 



■^'l.2,...,n+l 



(19) 



Again, one should contract all the indices corresponding to caps. 
3.1.5 General combinatorial formula for Kontsevich integral 

Now we know the contributions of all parts of the braid closure. In order to write down the answer, let us introduce the 
symbol Y\ representing the ordered product. 
The answer is then as follows. 

Let knot K be represented as the closure of a braid b E Bn: 



then the KI for the knot is given by the following expressiot^ 




(20) 



where the tensors T„, Xi and Bn are given by formulas (18), (13) and (19 1 correspondingly. 



3.2 Technique of computation 

In order to finally compute the coefficients of Kontsevich integral, one needs to take the following steps. 



First of all, note that basic elements of formula (20) are contractions T°-T'^ where one T" stands in some tensor component 



and the other in another one. Note that one can represent it as chord drawn on the knot from the place on the knot 
corresponding to the first to the place corresponding to the second T". Then, in order to compute the given order of 



(20), one has to consider all of the terms in (20 1 of that order, then for every such term draw a chord corresponding to every 



contracted pair of generators and untie the knot obtaining a chord diagram. 



•^One should understand that the symbol "tr" in \20\ stands for the contraction of the tensors Xi, Bn and Tn corresponding to their position 
in the braid. 
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More precisely, let n be a given order in h which we want to compute. In this order, Kontsevich integral is a Unear 
combination of group factors Gp = tr ^T°"f (^)r'^p(^)...T'^p(^")^ corresponding to pairings p = ^(ii, ji)(«2, j2)---(*n, jn)^ • 

One may associate with every such group factor a chord diagram, which corresponds to the particular pairing in a 
straightforward way. For example, there is the following correspondance: 



Thus, Kontsevich integral in order 71 is a linear combination of chord diagrams with n chords. 

Then, to obtain Kontsevich 



In these terms, formal formula (20 1 may be expressed as the following collection of steps. 
Let there be associators ^1, 



, and R- matrices Ri, . . . ,Ri assigned to a given knot K 
integral in order n one has to take the sum over all ordered partitions of n into k + 1 parts: 

n = 01 H h 0fc + ri H h n. 



(21) 



where all and are nonnegative integers. 

The R-matrix part is easier, for every R-matrix Ri we just insert consequential chords (propagators) in the corresponding 
place on the knot, and take the coefficient l/r^!. 



The associator part is trickier, there we have several summands due to the form of formula ( 15 1 for the associator, and to 



the fact that we should include also braids to the left, as it is done in formula (17). For every therefore there is a whole 
sum of ways to insert chords in the knot, with corresponding coefficients. That means that in a given order there are 



several terms of this order in expression (15 1. This number then multiplies even further when one takes into account (17). 

Then, for given partition n = 0i + • • • + 0fc + ri + • • • + r;, we take the product of all described above parts corresponding 
to R-matrices and associators. After expanding the brackets, we obtain a sum of several ways to insert n chords in the knot 
with corresponding coefficients. 

The only thing left is to untie the knot into a circle thus obtaining a linear combination of chord diagrams with n chords. 
Algorithmically, it can be done in a very straightforward way just by assigning an ordinal number to each part of the knot 
which belongs to an associator or to an R-matrix. 

Then we sum over all partitions of n and obtain the coefficients of order n in Kontsevich integral. 

Refer to Appendix |A] for more detailed description of the computation of KI coefficients for the figure-eight knot, i.e. 
knot 4i. 



4 Prom Kontsevich integral to Vassiliev invariants 
4.1 Chern-Simons definition of Vassiliev invariants 

As discussed above, the perturbation theory for CS leads to the following h expansion of vev for Wilson loops ([s]): 

00 dim(H„) 

<WH{K)>=J2fi'' ^n,-G„,™ (22) 

n=0 n=l 

here, K is the knot represented by the Wilson loop, R is the representation living on the loop, Vfe^m-are some rational 
coefficients depending on the knot and Gfc,m are the group factors called the chord diagrams with n chords. The chord 
diagrams with n chords form the vector space Hn- The coefficients Vn,m are not knot invariants because the chord diagrams 
are not independent, and to express this sum through invariants we have to choose some basis in 77„. The dimensions of Hn 
are summarized in the table: 



n 


1 


2 


3 


4 


5 


6 


dim{Hn) 


1 


1 


1 


3 


4 


9 



In order to pass to Vassiliev invariants we have to choose some basis in the space of chord diagrams. We do it following [TTj . 
refer to that paper for details. The so-called trivalent diagrams are introduced in a way represented for orders two and three 
in Figure [Sj Group-theoretical rules for graphical representation of chords and trivalent diagrams are presented in Figure [9j 
For the general definition of trivalent diagrams refer to [11], see also [12^. 

Let us explain the definition of trivalent diagrams on the first relation from Figure [8} 
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+ 2 







-H 3 



-+ 2 




Figure 8: Relation between trivalent diagrams and chord diagrams up to order 3 




^TrfT" ...) 



Figure 9: Group-theoretical rules 



rparpbrparpb rparpbrparpb _^ rparpbrpbrpa rparpb^pa^pb rparpb ^rparpb rpbrpa\ 



In Figure 10 one can find a collection of trivalent diagrams that form the so-called canonical basis {Gij} of H„ up to 
order six. 




S5I S?,j S?j3 












Figure 10: Trivalent diagrams 



Using this basis we rewrite (22) through invariants: 



00 dim(//„) 

,m Qn,r: 

n=0 m— 1 



(24) 



Here Vij are the so called finite-type or Vassiliev invariants of knots. They depend only on the knot under consideration but 
not on the group and its representation. 
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Now let us introduce the primitive Vassiliev invariants. It is a well known fact that the expansion of logarithm of any 
correlator in any QFT contains only connected Feynman diagrams (for more details about this situation in the Chern-Simons 
perturbation theory see [17]). This fact immediately leads to the following summation of 



OO Mn 

< Wn{K) > - n n i^'^Km Qn, 

n— m— 1 



(25) 



where Q'^ are connected diagrams, are primitive Vassiliev invariants. Here Mn is dimension of the space of connected 
chord diagrams (or equivalently the space of primitive Vassiliev invariants). The dimensions of these spaces up to order 6 
are given in the following table: 

n 1 9 s 4 .c; fi 

(26) 



n 


1 


2 


3 


4 


5 


6 




1 


1 


1 


2 


3 


5 



The meaning of the expression ( 25 1 is that V^.j in ( 24 ) are not independent. In fact only those coefficients Vij are independent, 
for which the corresponding diagram Qij is connected. Comparing h expansion of (251 with (24) we, for example, find: 



V4,l = 1/2 Vl,! 

And finally: The Vassiliev invariants form a graded ring freely generated by primitive invariants. 

4.2 Numerical results for invariants up to order 6 and families of knots 

Now one can use computational technique from the previous Section [3] and results from subSection 4.1 to calculate the 
Kontsevich integral. In this way we calculate Vassiliev invariants up to level 6 inclusive for knots with number of self- 
interSections up to 14 inclusive. These results are available in [18^ Here we list just a few examples: 





Va.i 


Va.i 


V4,l 


V4,2 


V4,3 


Vs.i 


V5,2 


V5,3 


V5,4 


Ve.i 


V6,2 


V6,3 


V6,4 


V6,5 


V6,6 


V6,7 


V6,8 


V6,9 


3i 


4 


8 


8 


62 
3 


10 
3 


32 


176 

3 


32 
3 


8 


32 
3 


32 


248 
3 


40 
3 


5071 
30 


58 
15 


3062 
45 


17 
18 


271 
30 


4i 


4 





8 


34 
3 


14 

3 














32 
3 





136 

3 


56 
3 


1231 
30 


142 
15 


1742 

45 


79 
18 


271 
30 


5i 


12 


40 


72 


174 


26 


480 


2512 
3 


448 
3 


104 


288 


800 


2088 


312 


41151 
10 


2494 
15 


7634 
5 


43 
2 


1951 
10 


52 


8 


24 


32 


268 


44 


192 


368 


64 


56 


256 


288 


2144 


352 


22951 


28 


29764 


137 


1351 


3 


3 


3 


3 


3 


15 


5 


45 


9 


15 


6i 


8 


8 


32 


116 


52 


64 


304 


64 


24 


256 


32 


928 


416 


2791 


884 


10084 


343 


871 




3 


3 


3 


3 


3 


3 


3 


15 


15 


45 


9 


15 


62 


4 


8 


8 


34 


38 


32 


208 


64 


24 


32 


32 


136 


152 


2129 


662 


1862 


463 


751 






3 


3 


3 


3 


3 


3 


3 


30 


15 


45 


18 


30 


63 


4 





8 


14 
3 


14 
3 














32 
3 





56 
3 


56 
3 


511 

30 


418 
15 


1858 
45 


65 
18 


449 
30 



(27) 

There is also one more technique to compute these invariants. It was suggested by Alvarez and Labastida in ^20J. We use 
this technique to check our results. 

4.2.1 Willerton's fish and families of knots 

On 18 Vassiliev invariant up to order six 

V2.I, V3,i, V4,l, V4,2, V4,3, V5,i, Vs.z, Vg.g, V5,4, Ve,!, V6,2, V6,3, V^A, Vg.s, Vg.e, Vg.T, Vg.g, V6,9 (28) 

there are 6 relations 



V4,i ^l/2Vl^, 


(29) 


V5,i =V2,iV3a, 


(30) 


^6,1 =^V|i, 


(31) 


V6^2 =1/2 Vl,!, 


(32) 


^6,3 =V2,lV4,2, 


(33) 


^6,4 =V2,lV4_3. 


(34) 
(35) 
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Remaining 12 Vassiliev invariants are generally thought to be independent. 

However, if one plots Va.i against V2.1 for all knots of given number of crossings, one obtains filled region with form 
resembling fish, which was discovered by Willerton in [21] and is hence called Willerton's fish, see Figure 11 See Appendix 



[B] where we present some new even more peculiar figures representing relations between Vassiliev invariants. 

The boundaries of the fish can be fitted with cubic polynomial. This suggests, that there are some additional relations 
on Vassiliev invarants. Note also, that for prime knots, which are actually plotted in Figure [TT] coordinates are "quantized" 
- they are integer multiples of V's for the trefoil. For higher Vassiliev invariants in our chosen basis it is not the case. 



! 1 I I I I I I 

! ! ! I I I ! ! 

■ t t 1 ! t t ■ 

i ! ! ! ! ! 




I : 

I ! 

i I 

. ! 



Figure 11: Willerton's fish, i.e. V3.1 plotted against V2,i for knots with up to 14 crossings 



Thus, it is natural to search for relations on Vassiliev invariants up to 6 order, i.e. to search for vanishing polynomials in 
V of order up to 6. If such polynomials can be found it would mean that Vassiliev invariants are not independent. However, 
it turned out, that there are no such vanishing expressions. 

Nevertheless, the following expression vanishes for all knots up to 6 crossings and for quite big number of knots with 
more crossings, for example, for 92 out of 165 knots with 10 crossings: 



Fi = 1/384 



16V2,i 
15 



2V2,i + V|,i - 6V3,i - 6V2,iV3a - 3Vi, - - 6V2,iV4,2+ 

+ 4V4,3 + 6V2,lV4,3 + 6V5,2 " 6V5,3 " 6V5,4 + 6Ve..5 - 6V6,6 - - 6V6,8 



(36) 



This is very interesting, because the general polynomial of order 6 has 17 free coefficients, while polynomial (36) vanishes on 
92 knots among solely the knots with 10 crossings! 

Another numerical experimental fact is that the value of ( 36 1 is generally close to the value of writhe number (see Section 
[5] for the definition), and tends to coincide with it, but not for all knots. 



5 Temporal gauge 

In Section [2] we declare that from Chern-Simons theory in Lorentz gauge one can get integral representations for Vassiliev 
invariants ^ [7| . Despite beautiful the integral representation for Vassiliev invariants is too complicated to compute and 
it is difficult to use them for investigations. From the other side, the knots are simple combinatorial objects and it is 
not surprisingly that simpler combinatorial formulas for the invariants should exist. In this Section we are going to define 
combinatorial objects related to the knots called writhe numbers. It will be shown that writhe numbers can be derived 
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from the Wilson loop operator in temporal gauge. However it turns out that it is impossible to express Vassiliev invariants 
through such writhe numbers. Because of this reason we define colored writhe numbers in a combinatorial way by analogy 
with the combinatorial definition of writhe numbers, because temporal gauge consideration does not lead us to them. Then 
it is possible to express Vassiliev invariants through colored writhe number up to order 4 at least (this fact was found in 

mm)- 

5.1 Writhe numbers from temporal gauge 

Consider Chern-Simons action in the temporal gauge Aq = 0. Then the propagator takes the following form: 

(AS(a;),4(y)) =0, fi = 0,1,2, (37) 
{A-^{x), Atiy)} = \^e^-8'^H{x^ - y^)5{x2 - 2/2)sign(xo - yo), ^ = 0, 1, 2. (38) 

Now let us consider the vacuum expectation value of the Wilson loop operator ([s]): 

{W{K)) = Y1 /'^^i /'^^2- f dxr,{A'''{xi)A^-{x2)...A'''{x,,))tr{T'''T''\..T''-). (39) 



Taking into account the propagators (37), (38 1 consider term of vev with ti = 2 in details: 

dx^, / dy^{A{Xf,)A{x^))) = 1 / / dx^,dy^e''''S''''S{xi ~ yi)6{x2 ~ y2)sign{xQ ~ yo). (40) 



2 _ 

Let us parametrize the knot by a parameter t running from to 1, then we can rewrite the last integral in the following form: 
' dhdt2 (^^~^^)s (xi(ii) - y.ih)) S {X2{h) - y2{t2)) sign (xo(ti) - yo(t2)) (41) 



Jo 



dti dt2 dti dt2 



To perform the integration we need to solve the following equations: 

xi{ti) ~ yi{t2) = 

(42) 

a^2(ii)- 2/2(^2) =0 

The solutions of these equations are the self-interSection points of two-dimensional curve {xi{t), X2{t)) which is the projection 
of the knot c on the plane {xi,X2)- Let us denote by ti < ^3 the values of the parameter t in the intersection points, then 
the two-dimensional delta-function in the integral can be represented in the following form: 

S{x,{t,) y,{t2)) S{x2{h) ~ y2{t2)) = )^ ,dxidy2 dx2dy,, 



dti dt2 dti dt- 
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Substituting this expression into (42) and integrating over ti and t2 we arrive to the following simple expression: 

/ dx^, / dy^{A{x^,)A{x^))) ^^^k, (43) 

where the quantities e^. are the "sings" of the intersection points. They can take values ±1 and are defined in the following 
way: 



dxi ^^^^ dy2 i^k^ _ dx2 ^^fc^ dyi 



* *' ' sign(x„(t;)-»((5)) (44) 



Thus, Chern-Simons theory in temporal gauge naturally leads us to the following notion. 

Definition 5.1. Writhe number for oriented knot w{IC) is defined as the sum of the signs of all the crossings: 



Wi 



(/C):=^e(p), (45) 
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Figure 12: Crossing signs 




Figure 13: Trefoil 



where the crossing signs arc +1 or —1 as indicated in Figure [T2j 

Writhe numbers can be represented graphically, if one gives another equivalent definition. Let us choose the origin on 
the knot and the orientation. When we are going along the knot, we meet every crossing point twice. We enumerate all 
crossing points by increasing sequence of the natural numbers. Then every crossing point is defined by the pair of numbers 



(11,12), ii 7^ i2 (see, for example. Figure 13). So, we get the following different definition of the writhe number. 



Definition 5.2. Writhe number for oriented knot w{IC) is defined as a sum of the crossing signs: 



Wi 



Now we can depict the writhe number by a chord diagram, which is called the Gauss diagram: 



Wi 



(^) = X! '^^i* 




(46) 



(47) 



Diagram language is more comfortable and illustrative for higher writhe numbers, definitions of which we give below. 

It turns out that the second Vassiliev invariant cannot be expressed through the just defined writhe number. Hence, we 
modify the definition of the writhe number. First, we introduce two types of the crossing signs e. 



Definition 5.3. 



0, 



if the strand, which led to the point (?i?2), is over another strand 
otherwise 



(48) 



£1112; if the strand, which led to the point (^112), is under another strand 
0, otherwise 



(49) 



Thus based on the last definition we can define two types of colored writhe numbers. 
Definition 5.4. Colored writhe numbers for oriented knot are defined as the following sums 



wT 



wl° 



OU 
«2 



EC 

l<i2 

/ J lit 



(50) 
(51) 



n<»2 
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Colored writhe numbers can be denoted by the following diagrams: 




Remark 5.5. We always read chord diagrams clockwise from the top. It is equivalent to introducing an origin point O 
slightly to the right from the top. 

5.2 Higher writhe numbers 

Any definition from previous subSection can be generahzed in a straightforward way. However because the colored writhe 



numbers are more general objects than the ordinary ones, we generalize only the last definition (5.4) 
Definition 5.6. Higher colored writhe number of n-th order for oriented knot w{IC) is a sum 



(53) 



tl<---<l2n 



where cti, . . . , (t„ take values ou or no. 



In graphical approach indices mi,...,m2n correspond to the points on the circle, corresponds to a line between 

points TOfe and m; (see the following example). 



Example 5.7. For order 2 we have only two diagrams: 



il<i2<i3<ii 



E 
E 



*i<*2<i3<U 




(54) 
(55) 



5.3 Relations between higher writhe numbers 

Higher colored writhe numbers are not independent. There are linear relations between them. These relations complicate 
the study of Vassiliev invariants because often it is difficult to reduce answers and understand their combinatorial context. 
For illustrative purpose we list here relations for orders 2, 3 and 4. 

5.3.1 Order 2 

There is the only relation: 




(56) 
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5.3.3 Order 4 



There are 195 relations and we do not list all of them. We just note that in order 4 relations do not mix connected and 
disconnected diagrams and diagrams with different topology. 

5.4 Vassiliev invariants via higher writhe numbers 

With the help of higher colored writhe numbers one can give very easy and elegant combinatorial formulas for Vassiliev 
invariants. 



5.4.1 Numerical experimental data 



From our numerical results (mentioned in Section 4.2 1 we find some useful properties of Vassiliev invariants. All of the results 
in this subSection are purely experimental ones. 

• Vassiliev invariants are can always be represented by sums of colored writhes with symmetric chord diagrams with 
respect to change of all over-crossings to under-crossings and under-crossings to over-crossings. It follows from the fact 
that Vassiliev invariants do not change under reflection. We call such combinations reflective. Therefore it is convenient 
to introduce reflective combinations of writhes as follows: 





• Vassiliev invariants can always be represented by sums of writhes with connected chord diagrams. For example, the 
first diagram in the figure below is connected while the second is not: 




Vassiliev invariants are always equal to a sum of writhes whose chord diagrams are irreducible, which means that we 
can not stick together any two chords in the diagram. Here we implicitly assume that to the left of the top point of 
the diagram we have the basepoint. Chords separated by the basepoint cannot be stuck together. For example, the 
first diagram in the figure below is reducible and the second one is irreducible: 




This property actually means that Vassiliev invariants of order k do not mix with invariants of order less than k. 
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5.4.2 Order 2 



One confusing fact related to writhes is their dependence on the choice of the base point on the knot projection. To construct 
invariants we should use only invariant combination of writhes that do not depend on such a choice. Invariant combination 
of writhes with two chords are the following ones: 




The second Vassiliev invariant should be some linear combination of these four expressions. We find: 




Note that due to the presence of the basepoint this diagram is irreducible. 
5.4.3 Order 3 

Invariant combinations of irreducible and connected writhes are the following ones: 
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We found the following interesting property of the first five combinations: 

C^i - C2 - C3 - C4 - C5 (80) 

For the third Vassiliev invariant we get: 

V3,i = ^ Ci + Ce (81) 

Acknowledgements 

We are grateful to A.Anokhina, E.Krylov, I.Polyubin, A.Popolitov, V.Slepukhin, Ye.Zenkevich and all participants of 
D.Anosov seminar for helpful remarks. We are also grateful to D.Bar-Natan for publishing part of our numerical results on 
the site katlas.org. We are especially grateful to A.Morozov for stimulating discussions. 

Our work is partly supported by Ministry of Education and Science of the Russian Federation under contracts 14.740.11.0081 
(P.DB., A.Sl.) and 14.740.11.0347 (A.Sm.), by RFBR grants 10-02-00499 (P.DB.), 10-01-00536 (A.Sl., A.Sm.) and 09-01- 
93106-NCNILa (A.Sm.), by joint grants 09-02-91005-ANF (P.DB., A.SL), 10-02-92109- Yaf-a (P.DB.), 11-01-92612-Royal 
Society(P.DB., A.SL, A.Sm.), by NWO grant 613.001.021 (P.DB.) and by Dynasty Foundation. 



A KI combinatorially for figure-eight knot 



In this appendix we will discuss the combinatorial technique of computing KI coefficients in more detail for the case of the 
figure-eight knot. 

Let K be the figure-eight knot, i.e. knot 4i. Braid representation for K with associators in correct positions is given 
in Figure 14 1. There are 8 associators 5'i,...,\E'8 and 4 R-matrices Ri, . . . , R4. Corresponding Sections of the knot are 
enclosed in thin dotted and dashed boxes. 

Let us describe the computation of KI coefficients of order 2. 

Recall that first of all one has to choose an integer partition of n = 2 into 12 parts (corresponding to R-matrices and 
associators). Note that there is no linear in T°T° term in ( 15 ), so all partitions n = (j>i + ■ ■ ■ + (fig + ri + ■ ■ ■ + r4 with any of 
4>i equal to 1 give vanishing contributions. 

Then we will have 4 partitions where both chords are taken from one of the R-matrices. One of the cases is drawn in 
Figure 141, where the propagators, i.e. chords, are represented by red dashed lines (the thicker ones). Further on, we will 



have 6 partitions where the chords are taken from two different R-matrices (Figure [l4p ) and 8 partitions where both chords 
are taken from one of the associators. This latter case is more complicated than the former ones, where we had only one 
term of a given order, according to (14 1. In the case of associators, however, we will have two complications: first, formula 



( 15 ) has already several terms of a given order and, second, the number of terms is further multiplied in accordance with 
(17) if we have lines to the left, as in the cases of ^'2 and '^r- Actually, for the order 2 we have just two relevant terms in 



(15), the one drawn in Figure 143 and the one which differs from it by interchanging their endpoints on the central string. 



Formula (17 1 implies that one should also include the same terms but with the leftmost endpoint shifted to the leftmost 



string, like in Figure 14 



The terms in KI corresponding to the choice of chords as in Figure 14 are as follows: 



1. Figure 

2. Figure 

3. Figure 
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r 24 ^ 
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Figure 14: Propagators 
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4. Figure fl4|3: {r^T^T^T^ 
C(3) 



5. Figure 



14 



'I927r3 



^^0,i rj-ia2 '^CLQ rpaQ rpaj rpa2 rpa \ rpaj rpa^ rpa4 ^ 



The final answer in order 2 for knot K is 



KI2 {K4J = ^ tr (t<^T''T''T'') - ^ tr {t''T''T''T^) 



(82) 



In orders higher than 2 the mixed terms where some chords are taken from associators and some from R-matrices start 



to appear. An example of such a chord configuration for order 7 is given in Figure 14 



B More pictures on Vassiliev invariants 



Here we list several more pictures of the type of Figure [TT] representing "relations" beween Vassiliev invariants. Very 
interesting is the form of the plot of Ve.s against V5,2, which is depicted in Figure 15 1. 



The following expression characterizes with good accuracy the long "horns" part of the figure by vanishing on most of 
the corresponding knots: 



-10752 



48616V2,i 
15 



178V2^i + ^2,1 - I254V34 + 66V2,iV34 - 3V| 1 

- 356V4,3 + 6V2,iV4^3 - 66V5^2 + 66V5,3 - 



- 356V4_2 — 6V2,lV4.2 — 
66V5.4 + 6V6.5 - 6V6,6 - 6V6,7 - 6V6,8, (83) 



see Figure [15) 3, where only points on which the expression vanishes are left. 

The plot of V^ 2 against V4_2 has vaguely the same form, and moreover the same expression ([83]) vanishes on most of the 
points in the "horns" part, see Figures [16^, [TBJa. 

Plots of 5V4_2 — 31V4,3 against V3.1 and V5,3 are also rather interesting, see Figures [17^, [TTja. 





Figure 15: a: horns, i.e. V6,5 plotted against V5,2 for knots with up to 14 crossings; b: same with constraint (831 
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a 



b 



Figure 16: a: horns, i.e. Vc, 2 plotted against V4_2 for knots with up to 14 crossings; b: same with constraint (831 



MM 



4 



Ji 




Figure 17: Birds, i.e. 5V4^2 — 31V4^3 plotted against Va^i (Figure a) and 5V4_2 — 31V4_3 plotted against Vs^a (Figure b), for 
knots with up to 14 crossings 
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